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The Matrix Product method (MPM) has been used in the past to generate varia- 
tional ansatzs of the ground state (GS) of spin chains and ladders. In this paper 
we apply the MPM to study the GS of conjugated polymers in the valence bond 
basis, exploiting the charge and spin conservation as well as the electron-hole and 
spin-parity symmetries. We employ the U — V — 6 Hamiltonian which is a simplified 
version of the PPP Hamiltonian. For several coupling constants U and V and dimer- 
izations 5 we compute the GS energy per monomer which agrees within a 2% — 4% 
accuracy with the DMRG results. We also show the evolution of the MP-variational 

parameters in the weak and strong dimerization regimes. 

PACS numbers: 71.35,74.70 



I. INTRODUCTION 

The study of conjugated polymers has been a subject of great interest for over two 
decades. There are both theoretical and technological reasons for this interest On 
the theoretical side, there exists a controversy within the scientific community over how to 
explain, understand and describe the photophysics/photochemistry of this class of materials. 
This controversy is of such a fundamental nature that the solution of the problem might 
be in a unification of the semiconductor and metal physics with the molecular quantum 
chemistry. On the technological side, pi-conjugated polymers behave as semiconductors and 
this has prompted several research groups to investigate the physics of these materials in an 



effort to determine their potential for improving the performance and efficiency and reducing 
the cost of hght-emitting diodes (LEDs). More recently, they are also considered to make 
an entrance in the field of photovohaics, where they could be used as solar cells. 

Saturated polymers are long chains of molecules, generally made of carbon with hydrogen 
on the sides, all attached to one another by single bonds. This constitutes the backbone 
of the macromolecule. The most relevant feature of these structures is the fact that the 
bonds are all single bonds or, in other words, that all the binding are of a- type. Saturated 
polymers are then all very insulating; they are not electronically interesting but are known 
for their flexibility although they are also quite mechanically strong materials. The most 
familiar of these compounds is the polyethylene. 

On the contrary, conjugated polymers show very interesting electronic properties together 
with remarkable mechanical properties; for instance, they can emit light and conduct elec- 
tricity In these compounds, two of the three 2p orbitals on each carbon atom hybridize 
with the 2s orbital to form three sp"^ molecular orbitals. These orbitals are responsible for 
the backbone of the molecular chain; these are the so-called a-orbitals. The third carbon 
orbital is Pz and points perpendicular to the chain. There exist a strong overlapping between 
nearest-neighbours Pz orbitals so that the corresponding electrons are fully delocalized on 
the whole molecule; these are the vr electrons responsible for all the interesting electronic 
properties of low-energy. For instance, because of these electrons, the linear chain, the poly- 
acetylene, which is considered in this work, is dimerized: its backbone shows an alternance 
between double and single bonds. Quite generally , despite a huge amount of works, the 
electronic properties of these compound stay rather controversial 0]. 

The delocalization character of the electrons in the vr molecular orbitals of the conjugate 
polymer chains led to the introduction of model Hamiltonians to study and predict their 
electronic properties. The initially simplest possible model is a tight-binding approximation 
or Hiickel model to describe the motion of tt electrons in a free way. This is a very 
crude approximation which has been improved in several fashions. One of them is the 
inclusion of electron-phonon interactions . However, this is not enough as the electronic 
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properties of vr-conjugated polymers derive from a true many-body problem were electron- 
electron interactions are equally important as the electron-phonon interactions. Then, the 
PPP (Pariser-Parr-Pople) Hamiltonian is used to model these electron-electron effects 

in a first approximation without taking into account phonon effects, to make simpler a first 
analysis of the electronic properties. In the PPP Hamiltonian, the alternating single-double 
bonds of the backbone polymer structure is realized by means of a dimerization term in the 
hopping kinetic energy. The most general form of the PPP Hamiltonian reads as follows: 

Hppp = Hk + Hi, (1) 
i/K = - E ^(1 - K-miclc,,^ + h-c.) (2) 
Hi = Uj:{clc,,-^){clc,,-'f) 

i 

+ E ^i^id) (E cU^i,^ - nei) (E CjVci,r - ^ol) (3) 

ij a- T 

where Hk is the dimerized tight binding kinetic part and Hi represents the Coulomb in- 
teractions among the electrons. Here the operators q^o-, cJ^ are standard creation and 
annihilation operators for vr electrons at carbon site i with spin a. The parameter t is the 
hopping overlapping integral between the nearest- neighbour carbon atoms, 6 measures the 
dimerization of the chain, nei, is the average number of electrons per site, U is the on-site 
Coulomb repulsion between the electrons, Vij is the distance between sites i and j along the 
chain and V{rij) is the long range contribution of the Coulomb repulsion. 

The PPP Hamiltonian has been the subject of extensive studies using a great variety 
of techniques such as Hartree-Fock, CI calculations, small cluster exact diagonalization, 
Quantum Monte Carlo and so on and so forth [||]. Only recently it has become possible 
to apply a new numerical technique, the Density Matrix Renormalization Group (DMRG) 
0, which allows us to obtain highly accurate results both for small, intermediate and large 



polymer chains [^12]. These DMRG studies have helped to clarify the correct ordering of 
excited states in the low energy part of the spectrum which are relevant for the nonlinear 
spectroscopic experiments. 
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In this paper we will concentrate on the study of the PPP Hamiltonian and leave the 
effect of interaction with phonons for future studies. The PPP Hamiltonian has been studied 
using an excitonic method based on a local description of the polymer in terms of monomers 
||13|| . The relevant electronic configurations are built on a small number of pertinent local 
excitations. This has provided a simple and microscopic physical approximate picture of 
the model. Recently we have extended this local configuration studies using the Recur- 



rent Variational Approach (RVA) method |Tj] in order to study larger polymer chains in a 



systematic way while retaining the previous intuitive physical picture. The RVA is a 



non-perturbative variational method in which one retains a single state as the best candidate 
for the ground state of the system. This reduction of degrees of freedom is initially done 
in order to keep the method manageable analytically. The aim of this analytical approach 
is to try to understand the relevant physical degrees of freedom so that we can figure out 
what the underlying physics is in a strongly correlated system. This initial analytical goal 
has been also developed in order to later acquire more numerical precision. To do this, the 
method becomes more numerical and somehow stands in between an analytical formulation 
of the DMRG and a numerical one. 

This effort of understanding the relevant electronic configurations in conjugated polymers 
has been also carried out in exact small cluster calculations using excitonic Valence Bond ba- 
sis 



p!q] for polymer chains of length up to 10 sites, arranged into diatomic ethylene molecules. 
A first comparison of RVA results with DMRG gave us promising perspectives to improve 
these variational calculations by incorporating more local configurations and variational 



parameters. In this paper we undertake this project by using a Matrix Product ansatz for 
the ground state (GS) wave functions This ansatz is a variational approach based on 
first order Recursion Relations (RR's) instead of second order RR's as in the RVA [|I8|. With 
this RR's we construct the GS of the polymer chain in different symmetry sectors based on 
the 16 local configurations of the diatomic ethylene molecule within the PPP approximation. 
Thus, the chain is built up by adding one ethylene at each step of the variational process. 
This paper is organized as follows. In Section 2 we introduce a Matrix Product ansatz 



specially adapted for the PPP Hamiltonian. In Section 3 we set up the Recurrent Relations 
to compute the GS energies in several sectors according to prescribed symmetries. In Section 
4 we present variational and DMRG results and make a comparison obtaining a very good 
agreement between them. Section 5 is devoted to prospects and conclusions. 



II. THE MATRIX PRODUCT ANSATZ 

The main idea of the MP method is to generate the ground state (GS) of a quasi-one 
dimensional system in terms of a set of states \a)N- generated by the following recursion 



formula 17,18 



I")a' = ^aAl^] \^)n \P)n-1 (4) 

where N denotes the number of lattice sites and Im)^- is a set of states located at the site 
N. For conjugated polymers each lattice site in refers to a monomer unit, and hence 
\m)N describes the 16 possible states associated to a single monomer. In table 1 we show 
the basis of local monomer states \m) n used in our construction. We have adopted a valence 
bond basis which is more convenient to our purposes although it can be easily related to the 
exciton-valence bond basis of references []T4| , |l6 . 



The states |a)Ar have to be regarded as block states made of intricate combinations of 
monomeric states whose structure depends of the MP amplitudes ^^^/^[m], which in fact are 
the variational parameters of the method. The latter parameters can be made to depend 
on the step of the RR, but in the thermodynamic limit one expect them to reach a fixed 
point value. Below we shall assume the thermodynamic limit, i.e. independence of y4Q,^^[m] 
on N, although computations can be done for any finite value of N. The choice of the block 
states I a) AT is mainly dictated by physical considerations and they are characterized by a set 
of quantum numbers as spin, charge, etc. In the case of conjugated polymers we shall keep 
6 block states which are to be thought as the GS's in the following sectors of the Hilbert 
space: i) singlet state at half filling with symmetry ^A'^, ii) singlet state at half filling with 
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symmetry ^B~, iii) a spin 1/2 doublet corresponding to making a hole to the half filled GS 
and iv) a spin 1/2 doublet corresponding to the addition of one electron to the half filled 
GS. The last two cases iii) and iv) describe localized charge transfer excitations between 
monomers, which play an important role in the GS of the polymer. In table 2 we give the 
6 blocks used in the MP-ansatz. 

Altogether we have a total of6xl6x6 = 576 possible MP amplitudes, but further 
constraints greatly reduce this number. First of all and without loose of generality one can 
impose that the block states |a)iv are orthonormal . This is guaranteed, for any value of N, 
by the following normalization conditions on the A's, 

AaA^] Aa',p[m] = 6a,a' (5) 

Moreover, the RR should preserve the charge and the spin of the states, refiected in 
the equations, 

ha = hm + hp (6) 
C2 _ cz _l_ qz 

where h^, hm, hp denote the number of holes and S*^, S^, denote the third component of 
the spin of the corresponding states. 

Finally, we can impose the conservation of the electron-hole and spin-parity symmetries 
generated by the operators J and P, whose action on a i*^- monomer is given by 0, 

J,|o) = -|x), P,|o) = -|o) 

J,|x) = |o), A|x) = |x) (7) 
J,| T) = (-1)^+^1 T), P^\ T) = -I i) 
U i) = (-1)''+^! i), Al i) = -I T) 

The action of J and P for a polymer with units is simply the tensor product of their 
actions on each monomer. In the eqs.(0) we use the convention according to which a state 
with symmetry ^A'^ has J = P = 1, a state with symmetry has J = — P = —1 while a 
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state with symmetry ^5+ has J = —P = 1 ( this differs in an overall sign to that used in 
0). The labels A and B refer to the reflection symmetry of the polymer, which shall not 
be imposed explicitely. 

Both the monomer states |m) and the block states |a) transform as follows under charge- 
transfer and spin-parity, 

J\a) = rii J\m) = rji^ \mj) (8) 

P\a) = 7]^ lap), P\m) = rj^ \mp) 

where 77;^ and r]^ can be derived from eqs.(|^, while t]^ and r]^ are the appropriated ones 
corresponding to the type of block chosen. In eq.(||) a j and mj denote the states obtained 
after the application of J on the states a and m respectively. All these quantities are given in 
tables 1 and 2. The MP equation (^) preserves the electron-hole and spin-parity symmetries 
provided the MP-amplitudes y4o,/3[m] satisfy the following constraints, 

Aaj,(3j [mj] = riir]ir]j [m] (9) 

Imposing the spin and charge conservation (^), the electron- hole and the spin-parity 
symmetries (|^) we are left with a total of 62 non vanishing MP-amplitudes y4Q,^/3[m] out 
of 576 possible ones. Moreover only 20 of these 62 parameters are independent. In table 
3 we give a choice for these parameters in terms of the MP-amplitudes, which we shall 
call hereafter Xi{i = 1, . . . ,20). Finally, the normalization conditions (|^) yield three more 
conditions on the set Xi given by, 

xl + xl + xl + 4x| + 4x1 = 1 

xl + x'^j + xl + 4x1 + 4a:?o = 1 (10) 
y-20 2 _ 1 

Hence altogether we are left with 17 independent variational parameters = 1, . . . , 17) 
which will be determined by minimization of the GS energy. Before we do that it is con- 
venient to parametrized the Xi parameters in terms of the yj ones (see below). From 
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physical reasons we expect that the most important MP-amphtudes will be given by 
Xi = y4ii[l],X8 = ^2,1 [3] and 0:17 = [9]. Indeed, Xi,Xs and xu correspond to the 
addition of a singlet, a local state, and a bonding spin 1/2 state to the GS block |1) , 
yielding a block state with the same type of symmetry as the monomeric state added. From 
this observation the parametrization we are looking for is given by 

xi = si, X2 = yisi, X3 = y2Si,X4 = ys^i, X5 = 

= y5S2, Xj = yQS2, Xs = S2, Xg = I/7S2, Xiq = ysS2 (H) 
Xn = y9S3,Xi2 = yioS3,Xi3 = ynS3,Xi4 = yi2S3,Xi5 = I/13S3 
Xie = yi4S3,Xi7 = S3,Xi8 = I/15S3, a^ig = ywS3,X20 = ?/l7S3 



Si = 1/Jl + yl + yi + 4yl + Ay: 



S2 = 1/Jl + yl + yi + Ay^r + 



s3 = i/^jE]Uy] 

The normalization conditions ( |10|) are automatically satisfied by the parametrization 
(|ll|), which on the other hand is quite convenient for numerical purposes ||T8[| . 

If we choose yj = 0(Vj) then the state |l)Ar generated by @ consists in the coherent 
superposition of singlets bonds on each monomer. On the other hand the RR's (^) also 



contain the Simpson state ||T^, which is the coherent superposition 



N 

ISimpson)^ = n (a;i|l)n + X2\2)n) (12) 

n=l 

With this state, the dimerized chain is viewed as a simple one-dimensional crystal of ethylene 
where, moreover, the electron correlations are ignored; this state was the reference state in 



1^ . It corresponds to yi = X2/X1 7^ and yj = O(for j > 2) 



III. GROUND STATE ENERGY 



In this section we shall briefly present the method for finding the GS energy of the MP 
ansatz whose minimization determines the MP parameters ( see references for more 
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details on the method). 

Conjugated polymers are customarily described by the Pariser-Parr-Pople (PPP) Hamil- 
tonian, however in our study we shall use a simplified version of it given by the U-V Hamil- 
tonian defined as, 

H = -tEi.Jl + + h.c.) (13) 

+ UY.i "■i.T'^ia + - ^)ini+i - 1) 

where and Ci^g are fermionic creation and destruction operators at the site i and spin 
s, Ui^s = c\gCi^s and = rij^-f + rij^j^. We shall work in units where the hopping amplitude 
t is set equal to one. The important parameters are therefore the dimerization 5, the on- 
site Hubbard coupling U and the nearest neighbour Coulomb interaction V . Since we are 
working in the monomer basis it is convenient to write the Hamiltonian (|l^) as follows, 

Hn = Ef=i h!f + Ef=-/ /^g+i (14) 

hf^ = + S]{clj_i^,C2j,s + h.c.) + U{n2j-i,^n2j-i,i + Ti2j,T^2i,i) 

+ \/K_i-l)K--l) 

^S-l = - S]{c2j,sC2j+l,s + h.c.) + V{n2j - l)(^2i+l " 1) 

where h^^^ is the intramonomer Hamiltonian of the j*^ monomer and ^^-^j+i is the inter- 
monomer Hamiltonian coupling the monomers j and J + 1. denotes the total number of 
monomers. 

The block states \a) n belong to different Hilbert spaces of the Hamiltonian (0), therefore 
the vacuum expectation value of will be diagonal with entries, 

=N {a\HM\a)N (15) 
The RR (^) yields a RR for these energies given by [0 



j^iN^i) = E + h^^^ + h^"^ (16) 



where 



To.,p = Em(^a,/3[m])2 (17) 

ha"* = 'Em,m' ,p AaAM^^A'^'] ei(m,m') 

(2) 

h» = Em's,/3,/3',7^a,/3["il]^a,/3'["^'l]^/3,7["^2]-4/3',7[m2] e2{mi, 7712] m[, m'^) 

ei(m,m') = (m|/i(i)|m') (18) 
€2(1711, m2;m[,m2) = {m2,mi\h^'^'>\m[,m2) 

The last two expressions are the intramonomer, i.e. ei, and intermonomer, i.e. €2, matrix 
elements in the monomer basis, which can be computed either analytically or numerically. 
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For the case of the PPP Hamiltonian, the number of these energy matrix elements (18 
huge and it is very lengthy the analytical computations of so many quantities. Instead, we 
have used numerical exact diagonalization techniques in order to compute them numerically 
once the PPP coupling constants are specified. This numerical coding is divided into two 
parts: 1) We construct the Hilbert space of states for the one- and two- monomer basis. This 
is done in a binary notation using a string of bits of length 4 for the one-monomers and 8 for 
the two-monomers. In the first half of each string of bits we encode the spin-up states and in 
the second half we encode the spin-downs. This representation we call it the tensorial basis. 
2) We represent numerically the action of the PPP Hamiltonian in the tensorial basis. This 
facilitates the computation of the energy matrix elements (^). Lastly, we perform several 
change of basis to bring the previous matrix elements to the Valence Bond basis employed 
in the variational recurrence relations. 



The RR (|T6D can be iterated to give once E^^ is known. Actually, the same is true for 



eq.(^) which gives the MP states |a)Ar once \a)i is given. We shall choose as initial states 
\a)i the lowest states of the monomer hamiltonian h^^^ in the corresponding Hilbert space 
sector. Hence the computation of E^ requires the diagonalization of ei(m, m'). 

Now the procedure goes as follows. Using eq.(|T6|) we find the value of E^ for a given 
set of variational parameters yj and look for the lowest possible value. This determines the 
value of these parameters and correspondingly that of the MP-amplitudes. One also finds 
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in this way the value of the GS energy density per monomer in the thermodynamic hmit, 

Coo = hm E^/N (19) 

IV. RESULTS 

In figure 1 we present the GS energies per monomer obtained with the MP method 
outhned above and the DMRG for the cases i) = 4, V = 1, 25 = 0.1,0.3,0.5, 1.5 and 
ii) t/ = 3, V = 1.2, 25 = 0.1,0.3,0.5, 1.5. For small dimerizations the relative error of the 
MP results as compared with the DMRG is around 4%, while for strong dimerization it is 
around 2%. 

In figure 2 we plot the absolute value of the 20 amplitudes Xi described in table 3 for weak 
dimerization {5 = 0.05) and strong dimerization {5 = 0.75) and couplings U = 3,V = 1.2 
in both cases. It is clear from fig. 2 that for strong dimerization the MP state is very 
close to the Simpson state for the most important amplitudes are Xi,X2,Xs and X17. For 
weak dimerization we observe a transfer of weight from these parameters to the remaining 
ones which show that the charge transfer excitations begin to play a more important role. 
This is specially clear in the behaviour of X4 which involves the monomer configurations 
(o I + I o), (o I + I o), (x t - T x)7(x i - i x), which are the typical local CT 
configurations appearing in the GS. On the contrary the parameter x^ remains very small 
showing that the monomer configurations (oT~T°))(°i~i°)5(x T + T x),(x | + | x) 
are very unlikely in the GS. 

These results are encouraging since they show that the MP approach gives a reasonable 
representation of the GS of the conjugated polymers in terms of a small number of variational 
parameters. They also show the possible improvements which can be achieved by first 
rejecting those monomer configurations which have small weight in the GS. One could also 
include blocks with spin 1 and singlet blocks with degeneracy. The latter type of blocks is 
needed in order to discuss the interesting crossing between the energy levels l^B^ and 2^A^ 
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V. CONCLUSIONS 



This paper represents the first attempt to generate a MP ansatz of the GS of conjugated 
polymers. Our results are rather encouraging since they show that we can get new insights 
and good numerical accuracy by improving the ansatz. Unlike other variational methods 
the MPM allows for a systematic improvement, becoming eventually exact when keeping a 
sufficient number of block states. Of course in the latter case the method becomes equivalent 
to the DMRG one [3T|. The usefulness of the MPM thus lies in a certain compromise between 
the desired numerical accuracy and the physical insight usually associated with the analytic 
nature of the method. The MPM also demands much less computing effort, an aspect which 
is certainly non negligible 
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FIG. 1. Ground state energies per site using the Matrix Product ansatz (blank square) com- 



pared with DMRG results (solid crosses) plotted against the degree of dimerization of the polymer 



chain and for several values of the PPP parameters U and V in equation (13). 
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FIG. 2. Absolute value of the twenty independent variational parameters used in the Matrix 

Product ansatz using six block bulk states. 
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Table 1.- States forming the monomer basis of the MP ansatz. The states given in the first 
column have to be normahzed. • — • represents a singlet valence bond state, x represents 
a double occupied site, o symbolizes an empty site and t, | symbolizes singly occupied sites 
with spin up and down, h denotes the excess or defect of holes as compared to the half 
filling situation. 2S^ is twice the third component of the spin, mj and mp are the states 
obtained upon applying the operators J and P on the monomer state m defined in eqs. 
d^). T];^ and r]^ are the corresponding signs appearing in eqs.(|). 
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Table 2. The notations are as in table 1. The states appearing in the first column are for 
illustration purposes. They simply show the type of symmetry of the block state as 
compared with the monomer states defined in table 1. 
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Table 3. List of the variational parameters Xj in terms of the MP-amplitudes ^^^^[m]. The 
total of non vanishing amplitudes y4a,/3[m] is 62. The remaining 42 amplitudes can be 

computed using eqs.(P). 
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